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Abstract 



We show how to perform systematicahy improvable variational calculations in 
the 0{2N) Gross-Neveu model for generic N, in such a way that all infinities usu- 
ally plaguing such calculations arc accounted for in a way compatible with the 
perturbative renormalization group . The final point is a general framework for 
the calculation of non-perturbative quantities like condensates, masses etc. . ., in an 
asymptotically free field theory. For the Gross-Neveu model, the numerical results 
obtained from a "2-loop" variational calculation are in a very good agreement with 
exact quantities down to low values of N. 
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1 Introduction. 



It has been advocated for a long time (|T]|-1T3I) that the convergence of conventional per- 
turbation theory may be improved by a variational procedure in which the separation of 
the action into "free" and "interaction" parts is made to depend on some set of auxiliary 
parameters. The results obtained by expanding to finite order in this "redefined" pertur- 
bation series are optimal in regions of the space of auxiliary parameters where they are 
least sensitive to these parameters. This is reasonable: such regions are those expected to 
resemble the most to the complete answer, where there should simply be no dependence 
on the auxiliary parameters. This intuitive argument can be made more precise, and 
there has recently appeared strong evidence that this optimized perturbation theory may 
indeed lead to a rigorously convergent series of approximations even in strong coupling 
cases. In particular, the convergence of this variational-like procedure has been rigorously 
established in the case of zero and one dimensional "field theories" [Tl[ . 



Having these rigorous proofs, it is very tempting to try to apply these methods in 
renormalizable quantum field theory. The main obstacle which one should overcome in 
this case is the compatibility of the variational expansion with the usual renormalization 
program of the theory. Indeed there are very few field-theoretical cases where the method 
has been tested and gave non-trivial results. Most of these cases concern the 0^ theory, 
which is unfortunately free in four dimensions, and the effective potential of the large- A?" 
Gross-Neveu (GN) model |, 0, |, [H. 



In this paper, in the finite-iV GN model, we shall show how one can build variational 
expansions which are compatible with perturbative renormalization and give non-trivial 
results even for low values of A^. Our approach can be in principle extended to more 
general theories, like QCD, as a new non-perturbative approach to the calculation of 
(dynamical) chiral symmetry breaking parameters |]T2| . 



We shall first take a closer look at the large- [|13| limit and from there infer the 
procedure for the case of arbitrary A^. Taking the fermion mass as a variational parameter 
we can obtain exact results in the large- A^ limit of the 0{2N) GN model, but as we shall 
see, in the finite- A^ case it seems necessary to go to infinitely high perturbative order just 
to obtain the right renormalization group (RG) behavior. We will see how to overcome 
this problem and interpret it in order to automatically include the correct renormalization 
group properties, obtaining finite answers as the space-time dimension goes to 2. This 
turns out to be a relatively simple exercise in renormalization theory, taken from an 
unusual point of view. At the end, we obtain a variational estimate of the vacuum energy 
and the mass gap, as a function of Ajjg, whose accuracy increases with the order of the 



perturbative expansion. Exact results for these two quantities are known [[14| , [15| and 
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provide a test for the accuracy of the method. Note that our method also has definite 
connections with previous derivations of a mass gap in the framework of the (non-hnear) 
sigma model |TB|, In particular our way of using the RG properties to relate the 



mass gap to the basic scale, A;^, is quite analogous with the one in the latter reference, 
although the subsequent method to extrapolate the pure RG results to the actual mass gap 
values is completely different. Moreover, as above mentioned the essential novelty in the 
present approach is to clarify the link between RG results and the variational expansion, 
allowing a treatment of the latter that is consistent with infinities and therefore much 
more adaptable to other theories. 

In section ^ we explain the basic idea and ingredients of our method, using only the 
first coefficients of the renormalization group (3{g) and 7(5') functions in order to keep the 
discussion clear. In section ^, we show how to generalize, formally, to arbitrary higher 
orders the previous construction, and derive explicit formulas for the fermion mass and 
the vacuum energy density using the presently known maximal perturbative order, i.e. 
the three-loop renormalization group functions. In section ^, we construct an alternative 
Ansatz, starting directly from renormalized quantities, which is shown to differ in general 
from the previous one by a specific renormalization scheme change, and lead to essentially 
similar results. It also gives a more transparent construction, in particular more appro- 
priate to generalization to higher orders and other theories. Finally some technicalities 
are treated in two Appendices. 



2 One-loop renormalization group properties 



Let us start from the 0(2N) GN Lagrangian |T^, modified with a non-zero fermion mass 
as follows ^ 

^G7v(mo, go. x) = ^{i^^d^^ - mo)^ + x {^{^^)^ + mo^^); (2.1) 

where we introduced for convenience a parameter x, interpolating between the free, mas- 
sive Lagrangian for x = 0, and the usual massless GN interaction Lagrangian, for x = 1 
(summation over the 2A^ component of fermions is implicitely understood Q) • As is well 
known, one can solve the model exactly in the large N limit and obtain a non-trivial mass 
gap in the massless limit, i.e. for x — 1 in the present context. Now how can we treat 
the theory as defined by ( p.lD for arbitrary A^? A direct application of the arguments 
developed in the introduction would consist of doing perturbation theory with respect 

^We consider the 0{2N) case only for convenience of comparison with the resuhs of ref.||T^, but in 
fact there are a priori no hmitations in our method to consider N half-integer as well, i.e. corresponding 
to the usual 0{N) model. 
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to the reorganized Lagrangian above, i.e. expanding the above terms to some order in 
X. One would then look for extrema with respect to m, supposedly approximating the 
unknown exact result for a; — > 1, which by definition does not depend on m. Lessons from 
the anharmonic oscillator case leads one to expect a systematic improvement when going 
to higher orders in the expansion in the better approximation being assumed for the 
fiattest extrema structure. 

Unfortunately, one immediately encounters a number of obstacles. Most importantly, 
before accessing to any physical quantity of interest for such an optimization the theory 
has to be renormalized, and there is an unavoidable mismatch between the expansion in x 
as introduced above and the ordinary perturbative expansion as dictated by the necessary 
counterterms, as we shall see explicitely in section 2.2. Independently of that, even in 
the most optimistic case where arbitrary orders of perturbation theory would be known 
(which is only true in the large N case), it is easily seen that at any finite order one 
only recovers a trivial result for x ^ 1 (i.e. mo — >■ 0), which is the limit we are however 
interested in to identify a non trivial mass gap. 

Actually these facts can be circumvented by advocating a certain resummation Ansatz, 
whose main ingredient exploits the renormalization group invariance of the theory and 
analytic continuation properties of the (arbitrary) mass parameter m ^, which at the end 
defines a mass gap as an integral transform as we shall derive below. This (renormalization 
group invariant) Ansatz is exact in the large limit, where it can be shown to resum 
the X dependence exactly. In the finite N case it is only exact as far as the leading (or 
next-to-leading) renormalization group behavior is concerned. But the important point 
is that it also provides a non-trivial resummation of the perturbative serie in x, in the 
sense that it gives a finite and non-zero result in the m — limit. This is expected to 
be a sufficiently good "trial" Ansatz for a subsequent optimization with respect to m, as 
will be motivated below. 



2.1 The mass gap in the large N hmit 

To be more concrete let us start with the expression of the one-loop mass and coupling 
counterterms of the model, expressed with the renormalized parameters m and g related 
to mo and Qq by : 

m 

mo = mZ„ 



[1 



^More precisely, of an arbitrary parameter related to the original Lagrangian mass mo in a way to be 
specified next. 
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gl = g'i.-^Z, = ^^^T) (2-2) 

where fi is the arbitrary scale introduced by dimensional regularization. Note that those 
one- loop expressions are exact in the N ^ oo limit, where it can be shown that there 
are no other leading corrections. For finite this one-loop RG dependence generates 
however only the leading terms in 1/e, to all orders. 

Consider now the following expression for a (bare) resummed mass, as inspired by the 
latter RG properties: 

mpimo) = ^ ^TTTF (2-3) 

(1 - 2{N - l)g^o^fi-'mUmo))^^ 

(where T = r[— e/2]/(47r)^"'"^/^). Eq. (|2.3|) is obviously renormalization group invariant 
since expressed only in terms of the bare parameters mo and go, and can be expressed 
in terms of the renormalized parameters m(/i) and g{fi) with the help of eq. ( p.2| ) above. 
When this is done, one obtains for mp d.n expression which can be recursively expanded 
in powers of 5'^(/i) and, as easily checked, is finite at each order as e ^ 0, at fixed g'^{fi) 
and m(/i), thanks to the recursive dependence in m^? in eq. (|2.3| ). Explicitely one obtains 



mF(m) = m [1 + -f In : ^ ]-^^ ^ 2.4 

where the MS scheme scale, introduced for convenience, and 

g = g{fi), rh = m{jx). Eq. ( p.4|) resums the leading logarithmic dependence in m, and in 
the iV — > oo limit there are no other corrections. Using the recursivity one can rewrite 
( p.4|) identically as 

m^" 

m^ (m") = ^3T7F ; (2.5) 

f{m") 'v-i 

where the scale-invariant, dimensionless mass parameter: 

iV-l/2 

m" = [ni/Ajjg) [{N — l)g /tx] ^-^ has been introduced for convenience and is related 
to / as 

TV — 1/2 

f{m") + ^ _ / In /(m") = In m" (2.6) 



while ^yjig = [X exp[—T{/{{N — l)g'^{fi)) ] is the MS RG-invariant basic scale. 

N-l 

Now from (|2.6| ), it is easily seen that for m" — > (i.e. m — * 0), f{m") ^ (m")^-i/2^ so 
that one recovers the well known result, 

= Ams (2-7) 

for the mass-gap in the large A^ limit. 
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2.2 Usual plague of the x expansion and its cure 



Up to now we have not made any use of the expansion in the "variational" parameter x 
and indeed, as expected, that was not necessary to obtain the correct result for the mass 
gap in the N oo limit. Since the large N result is a good consistency check, and the 
perturbation in x will play a definite role later on in the more complicated arbitrary N 
case, let us see what is happening when considering the x expansion. 

Introducing the x expansion parameter in formally amounts to substitute ttlq 



mo(l — x), (7q — >■ QqX in any (bare) expression, in particular in (2^). One however soon 
realizes that ( p. 3D is then no longer finite for e ^ at any order, apart for x = 1, 
in which case only the trivial result mi;' = is recovered, as announced. Indeed the 
resummed expression (|2.4| ) resulted from the usual RG properties and has no reason to 
be compatible with the way the perturbative x expansion is introduced. The cure to 
that problem is to resum the series mp{x) = Y^'^=o^q^'^ i ^ different way: this is made 



possible by analytic continuation in x, and an adequate contour integral can be shown [13 



to resum exactly (in the large N limit) the above series, as explained in some details in 
Appendix A. Apart from technical details, the main effect of the resummation is to cancel 
the factor of {1 — x) which arises from the substitution applied to uiq and was responsible 
for the trivial result at any arbitrary order for a; — >■ 1. The net result gives the mass gap 
as a specific contour integral expression: 



mp{m'Q) 
with 



mo 

2iTT 



due^ 



iV-l/2 



2{N-l)glT {m'.uY [/i («)]"- 



iV-l/2 



(2.8) 



where u = q{l — x) has been introduced as an appropriate change of variable to analyse 
the X —>■ 1, q oo limit of the g*^-order expansion of mpix), ttIq is obtained by the 
rescaling uiq = q m'g, and the integration contour runs counterclockwise along the cut 
negative real axis. 



The other nice thing about equation ( |2.8| ) is, that it has a smooth limit as e ^ 0, at 
fixed renormalized m' and g'^, where m' is obviously related to m'o as in (|2.2| ). Indeed in 
the limit e — > 0, we have : 

mF{m' 

^^^^ [f[iu)]^ 
with: } (2.9) 



m 
2i7r 



m'u 



n 



A^- 1/2 
A^- 1 
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As expected, this can be expressed in terms of the dimensionless parameter m", given by 



m' 



m" = ^ ^ , (2.10) 

which is invariant under the renormahzation group of the massive theory i.e. 

r d g\N-l) d g\N - 1/2) , d . 

g- m-—\m =0 . (2.11) 

^ Ojj, TT og TT am' ' 

After these rescahngs one obtains: 

mF{m") = Aj^ — f , (2.12) 

withQ 

/(n) = ln(m%)-^^^ln/(«) . (2.13) 

This conceptually very simple expression defines a function mFlm"), which should be 
studied, looking for extrema, and/or its value at m" = 0, as explained in references |I3|,P^. 



In fact, as far as the large N limit is concerned the optimization is trivial, as one may 
expect: from (|2.13| ) one can find the series expansion of f{u) near the origin as f{u) ~ 



{m"u) ^ from which it is immediate that the expression ( ^^.12] ) has a simple pole at 
M = 0, whose residu simply gives mp = A-j^. Now however, the main interest of the above 
construction is that the resummation of the x expansion remains valid, when considering 
the generalization at higher RG orders of the different expressions, and when one includes 
as well the non-RG perturbative corrections which are present at the next order, as we 
shall see. Remark at this point that all the complexities of the renormahzation procedure 
are hidden in the definition of the function f{u) in ( |2.13| ). We note that once the pure 
number m" is fixed, for example at an extremum of mp{m"), the corresponding value of 
miT' is an invariant of the renormahzation group of the massless theory, as is evident from 
equation ( |2.12| ). 



3 Higher order results 
3.1 Generalized Ansatz 

For finite A^, and at a given order in ^f^, one knows Zm and Zg from the knowledge of the 
renormahzation group functions j3{g) and 7(5') at the appropriate order in g^. Hence we 
can ask the question : 

^Note that f{u) is related to f[{u) in (U) as f{u) = -^^^^Ji^ f[{u) . 
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Can we guess for mp (or any other physical quantity) an expansion in powers of (y'o'^o 
which is finite when £ — at fixed and m and has a power expansion in g"^ ? 

The exact expression for mp is of course an answer, but certainly not the only one, to 
this simple purely mathematical question. What we are after is a class of answers, and a 
general principle which will enable us to pick in this class the appropriate answer which 
reproduces the results of the calculation up to a given order in the loop expansion, i.e. 
in the expansion of in powers of g^m^. The graphs of figure |l| give the perturbative 
expression for mpimo) to (9((7q), for arbitrary A^: 

mF(mo) = mo + (2Ar - 1) ^2 f 

+ (2iV - 1) gt m^^ (1 + .) { (2iV - 1) - } 

+0{gl) , (3.1) 

where we have used the following definitions: 

f E 



r(-|) 

(4vr) 1+^/2 



^ r(-£/2) k \{\-yyx(\-x)^y\ 



3£ 



1 + 

£^0 2 



2 2.9022... 



£^ + 0(£^) . (3.2) 



-24 2 

Now the previous expression (^73|) (from now on denoted as JF), 

^("^o) = ^TTTT (3.3) 

{\-2{N -\)glV^--T^}-^ 

reproduces eq. (|3.1|) only to order g^ and the leading, l/e^ divergent terms of order 
(^o; while at the same time being finite, to all orders in g^ fixed for e — > 0, when the 
renormalized mass and coupling are substituted, as discussed in the previous section. It 
should be noted however that, even when restricting to the first RG order, mp need not 
be identical to JF : at any order in g^ there is room for a finite term, to be determined by 
an explicit calculation. This follows from the fact that 

^^^1-^ (3.4) 

is itself, as well as all its powers, in the class of renormalization group invariant functions 
of g^ which are finite as e — > 0, at fixed g^ . So it is clear that mp can be expressed as a 
unique expansion in this expression. The above fact enables us to make compatible our 
variational procedure with renormalization theory at finite as we shall see next. We 
must accordingly take into account the fact that equation (|H73|) does not reproduce the 
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expansion of equation at order g^, since a finite term (as e 0) survives, as well as 
a next-to-leading divergent term ~ Qq/e, at that order 0. 

To systematically correct these discrepancies, while keeping at all times smooth limits 
as e ^ for fixed renormalized parameters, is achieved by noting that the renormalization 
group invariant: 

~ ~ , , iv-i/2 (3-5) 

is finite for e —>■ and has a power series expansion in both and g^. With this invariant 
and all its positive powers, one will then reproduce the whole perturbative expansion (in 
powers of gQ ) of equation (|3.1|). 



The construction of our generalized Ansatz therefore goes as follows. We consider 



entering the following renormalization group invariant expression for mp, 

mo 



rriFimo) 



Xg^ml 



P 90^1' ^ 



{i-2(iv-i)^?gr(^r} 



(3.7) 



where the constants A, B and D are to be determined by requiring equation (|3.7|) 
to be finite as e — > when expressed in terms of m and g^. A, B and D are therefore 
obtained as an expansion in powers of e, where the different orders of the RG coefficients 
enter at the different order of e. The coefficient X, p, ■ ■ ■ are then given by matching 
the perturbative expression for mp, equation (|3.1| ). This is explained in more details in 
Appendix B. 



With expression (p.7|) and a straightforward generalization of the steps decribed in 
Appendix A, one can reach infinite order in the variational expansion in a very similar 
way as the one leading to equation (p.9|), ( |2.12| ), and obtains the following generalized 
version of equation ( p.l2| ): 



MS ^ 



m" f due^ r Ml M 
77— f N-i/2 11 + -r + ~7 



* Indeed, even at order g^, a finite term survives, if in equation (2.3) we replace r(— |) for example 
by its leading behavior, — 2/e. 
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where now I 



f2{u) 

Ml 
M2 



\n{m"u) — 



N 



N 



3(A^- 1/2) 

{N - l/2)(36iV2 - Q2N + 17) 
64(Ar- 1)4 

-^-^^^-0.725551... + ^) 

2(A^-1)2^ 48^ 



(3.9) 



and in the overall factor in ( p.8|) we again identifed the usual renormalization group 
invariant scale at this order, 



A- 



MS ~ ^-^P 



vr 



{N - l)g^ 



271 



1 

2{N-1) 



1 



y \ 2{N-i) 



up to an overall factor independent of in (p.8|) . 



We now discuss expression (3^), which is the essential result of this paper. One can 
analyze the function mi?(m") in standard ways. One finds that for m" — > 00, its behavior 
is accurately described by the first term of equation 



A (-] 2{iv-i) 



m 



iv-1/2 ' 

(Inm") 



(3.10) 



a typical renormalization group result, and that the A^„//^ correction contributes a 
corrective term of order (Inm")"*^ to this, as well as In Inm" terms. Hence, by going to 
higher and higher orders of ordinary perturbation theory, one determines the function 
mi?(m") with an increasing accuracy for m" large enough. 

The same phenomenon had been noticed for the anharmonic oscillator case in ref- 
Oi- Just as in that case, the hope is then that already with a low order 



erences 
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calculation, one may reach an accurate value for, say, the physical mass gap by taking an 
extremum of mp{m"), the best one being presumably the one closest to m" = 0. For the 
N = 00 case, we have seen that all the coefficients are zero, and the extremum at 
m" = is the exact answer. For finite N, the numerical convergence of the approximation 
can only be decided by an explicit calculation. 

Notice that, when the order of the perturbative expansion increases, one can deduce 
from the behavior of f{u) as u goes to zero, that mF{Tn") becomes more and more singular 
as m" — > 0. This same fact had been noticed on the case of the anharmonic oscillator. 



'Note the difference between the coefficient of the ln/2 term in (3.9) and the power of /2 in (3.8). The 



previous equahty of those quantities, in eq.(2.£) (and the consequent simple pole behavior for m" — s- 0) 
was an accident of the first RG order. 
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There, this singular behavior at m" = did not prevent the analogous function from 
converging to the right answer for a wider and wider range of the variational parameter 
when the order of the expansion increases. Thus, one could obtain very accurate results, in 
particular when resummation techniques are used to improve the behavior of the integrand 
of equation ( p.8|) for m ^ 0. This will be discussed below. 



3.2 Vacuum energy density 

A similar treatment can be applied to the vacuum energy density , with an appro- 
priate modification to account for the structure of its divergences. 

Indeed, in contrast with rriF, the vacuum energy, in the massive theory, is not multi- 
plicatively renormalizable. Rather, is, and consequently is finite. Combining 

^ ( Fi) 

this fact with dimensional arguments, it follows that all infinities of Eq which persist, 
at a given order of the perturbative expansion (after coupling constant, mass, and wave 
function renormalization), are proportional to Hence, by subtracting from Eq^^ the 

quantity Y^h{e) with a suitably chosen function h{e), Eq^'' can be made finite for e ^ 
at fixed renormalized coupling constant g and mass m. It is crucial to remark that such 
a subtracted counterterm gives vanishing contribution to our procedure for any e as soon 
as the order of the variational-perturbative expansion is larger than one0. Expanding now 
the regular function h{e) as 

h{e) = Y.Ke^ , (3.11) 

n>0 

one can uniquely determine the coefficient of the term by a direct perturbative cal- 
culation to n + 1 loops. In contrast with the mass case, the fact that the {n + 1) loop 
information is needed to determine unambiguously the 1//" perturbative corrections in 
our Ansatz, is of course a reminiscence of the original ambiguity due to the above men- 
tioned additional divergences of the vacuum energy []. 

Once this is done, the finite terms that are left over can be treated in the same way 
as mp was treated. We skip the details and only give the final answer, as a function of 
the variational parameter m", at e = 

e 



Eo{m") = 2(iV- 1) A|^ (3.12) 

,//2 



m" f udue^ r E\ £• 



2m J fN-i ( \ f2 f: 



C 01 Oo 



^ More precisely the subtraction procedure does not give any new contributions beyond the generic 
form in ( |3.12 ) below, but actually the specific values of the perturbative coefficients, £i in ( 3.f 2 ), consis- 



tently include a dependence on the subtracted terms. 

''Actually however, only the value of the simple pole in e at order n + f is needed. 
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where 

Si 



N{N^ + N/2-?>/A) 
2{N - l/2)(iV- 1)2 



_ N{N-l/2){N'-QN + 7/2) 

~ 8(iV-2/3)(iV-l)4 ■ ^^-^^^ 



3.3 Refinements 



While our main result formulas (|3.8| ) and ( p. 121) contain in their derivation all the concep- 



tual ingredients of the method, they leave room for some refinements, which are necessary 
before they can be of practical use for optimization. 

The first one comes from realizing that, even at a given RG order, our different Ansatze 
are not unique, since there are infinitely many ways of introducing the purely perturbative 
terms, 1//" corrections. (Moreover beyond the first RG order even the resummation of 
the correct two-loop RG behaviour is not uniquely fixed by the function JF as defined in 
eq. (|0|), as will be discussed in the next section). To take into account this freedom, let 
us introduce an arbitrary parameter a, from substituting ft ^ a fx in the different relevant 
expressions. This constant, which enlarges the class of our renormalization group invariant 
finite guesses for mi;' , cannot be ruled out a priori, and parameterizes a renormalization 
scheme (RS) dependence. Of course, if one were able to work to all orders, the final 
answer would be independent on a. However, since we are truncating the expansion in 



equation (|3.7| ), the choice of this scaling constant a turns out to be important in order 



to obtain a reliable numerical estimation for mp. In the present model, the value of 



a which ensures a rapid convergence of the expansion in equation (|3.7|) is expected to 
be a = 1 + 0{1/N). We have however no a priori idea of the value to chose in more 
complicated theories as QCD. Hence we shall see in the following how to cope with the 
existence of the parameter a, intrinsic to renormalizable theories, in order to avoid too 
much additional arbitrariness in our results. 

Applying now the variational procedure, we finally obtain a last version of equa- 
tion (|3.8|) , expressed in terms of the usual renormalization group invariant parameter 

/e\'2uti) m" a r due"" c Mx{a) M2{a) ^ 
^Ms 2nr J jTv^(^) ^ J2 I2 ^ 

where the recursively defined function /2 is given by 

. ,3.15) 
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The coefficients in the expansion are now given by 

{N -1/2){N^ -5/2N + 9/8) 

(iV- 1)3 

where A^i(l) and A^2(l) are defined in equation (|3.9| ). 
Similarly, we obtain for the energy 

Eo fe^-^^ {m"fa? r udue"" r Si{a) ^2(0) , \ 1 7n 



In a , (3.16) 



with 



£^1(1) and £^2(1) being given in equation ( p.l3|) . 



We are now asked to find a stationary value for these expressions in term of m" : 

c)m" dm" 

In order to get a feeling of what we have achieved, we have plotted in figure ^ the 

values of the mass gap mp/Aj^ versus the variational parameter m", for different values 

of the scaling constant a ranging from a = 1.1 to a = 1.2. Similarly, we have plotted in 

figure § the values of the vacuum energy Eq/Ajj^"^ versus m", for a ranging from a = 0.9 

to a = 1.1. In both cases, we have taken a reasonably small number of "fiavors", = 3. 

We have limited the graph to the vicinity of the origin in m" since it is the region of 

interest (for large m" the curves are following the expected perturbative behavior). The 

dashed horizontals represent the known exact values for the mass ^ |]T^, obtained from 

the Bethe Ansatz: ^ 

niF (4e)2(N^ 



^MS r(l - 2(77ZT) 



(3.20) 



®In these figures, we have used an expansion for m" — > in order to produce the curves in a reasonable 
computing time. The corresponding sHght change in the shapes of the displayed curves is irrelevant for 
the discussion. 
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and for the vacuum energy[| 



^° = -X"*2(iv3T) • (3,21) 

Many comments are now in order: 

First of all, we have succeeded in producing finite values for the mass gap mp and 
the vacuum energy Eq at 2-loops order, including a variational mass m", in a way consis- 
tent with the renormalization group . The extension of the procedure to higher orders in 
perturbation is only a matter of computational effort. Moreover, the method can be ap- 
plied, at least in principle, to some of the order parameters governing the chiral symmetry 
breakdown in the realistic case of QCD, as we shall see in a subsequent paper. 



Secondly, these expressions, equations (|3.14|) and (|3.17|) , exhibit stationary values with 



respect to the parameter m", as can be seen on figures |^ and ^ . We can then proceed 
and look at these stationary values. 

We immediately meet an obvious problem: due to the large dispersion of the curves 
when the scaling constant a is varied, it seems difficult to predict a reliable numerical 
value for the ratios mp/ and Eq/ Aj^g. Actually, the stationary values for the ratio 
nip/A-jjg are seen to lie in the range 1 to 3 when a is varied from 1.1 to 1.2, and one 
cannot hope to extract from that more than a vague estimate of the mass gap. 

However, one can see in figure ^ that the flattest curve, whose stability plateau gets 
closest to m" = 0, lies quite near to the expected exact value. Moreover, the divergence 
of the curves near m" = 0, where we would expect to obtain the exact answer, as in the 
N ^ oo case, is only due to the breakdown of perturbation theory as m" — * 0. This, 
and previous experience with the anharmonic oscillator case [jl9|, suggest us to try an 



extrapolation towards m" = in order to improve the numerical precision of our results. 
This will be done now. 



3.4 Optimization results 



The strategy is to turn the expansion in powers of I//2 of equation ( |3.14 ) into an expansion 



around m" = 00, whose first term appears in equation ( p.lO|) , later to be extrapolated 



by some variant of Fade techniques to m" = 0. This however produces an expansion 
containing Inln. . .Inm" , which is not suitable for an extrapolation towards m" = 0. In 



'We are grateful to Prof. Al. B. Zamolodchikov for communicating his result to us before publication. 
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fact, equation ( |3.15|) suggests that instead of m" one should use the variable rj defined by 



V 



In terms of 1], the Inln . . . Inm" disentangle, and expression ( |3.14| ) now admits a very 
nice asymptotic expansion for rj ^ oo, which is readily seen to giveQ 



Am5 iV-1 ' ^ V iV-l/2 

Taking now the derivative of the logarithm of this expression yields a power series in 
l/t] , which is suitable for a Pade approximant analysis: 

dri 7] N -I r/2 V 1/2 ^ ; ^ 

(2iV-3/2) J^. 

TV - 1 iV - 1 ^ ^ ^ 

From the large-iV case studied in a preceding paper, we can make the guess that, 
were it not for the fact that perturbation theory breaks down at m" = 0, mF{m") would 
approach the exact mass gap as 

rriFim") ~ m^(0) + 0(e-^°°^*/"^") . (3.25) 
Therefore we can try a Pade approximant of the form 



dr] ~ 1 + C^h + Cs/r/' + . . . ' ^ ' ' 

which ensures an exponential convergence of mp near m" = 0, provided that Ci < 0. 
However this is by no way limitative, in particular for finite N there are no guarantees that 
the exponential behavior persists, therefore we shall also use simple Pade approximants 
(i.e similar to ( |3.2(j| ) without exponential factor), and compare different Pade results. For 
any Pade approximant type we can now determine the coefficients { Ci , C2 , C3 . . . } by 
matching the large-77 expansion (|3.24| ), and integrating equation (|3.26|) from rj +00, 



IOt 



In (3.23) the cocfBcicnts of r] " involve integrals of the type I{p) = § dy In^ y, which are easily 



evaluated as a (finite) expansion of real integrals expressable in terms of Euler Gamma functions. 
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where mp is given by equation ( p.24|) , to 77 = 0, will produce an extrapolated value at 
m" = for the mass gap mp/Ajjg |^. 

Very likeky, it should exist an optimal Pade order, which may however vary as 
varies, and is difficult to guess a priori. One may naively assume at ffist that the best 
Pade approximant should be of the same order of the known perturbative terms. Ac- 
tually the situation is slighlty more complicated, since the RG coefficients, which is a 
genuine information independent of the purely perturbative terms, also enter Pade ap- 
proximant expansions. Moreover, the larger is, the more perturbative information may 
be considered available, since all the purely perturbative coefficents (9(1//") in mp are 
vanishing for N —>■ 00 (although we do not know how fast these are vanishing as functions 
of N). Hence, one can test the convergence properties of the increasing order Pade ap- 
proximants ( p.26| ) in the large- A^ case, and infer from it the accuracy of the lower Pades 



at finite A^. Without much loose of generality in our numerical study we compare Pade 
approximants of different orders, from (1,1) to (2,3). 

The results of the procedure for N 00 are presented in table where we have given 
for different values of the scaling constant a, the value of the extrapolated ratio, using a 
Pade 1-1, a Pade 2-2 or a Pade 3-3 approximant^. The exact value is mp/Aj^ = 1. We 
see that the stationary (in a) values of the successive Pades do indeed converge, albeit 
rather slowly, as the order increases, and with what may be considered relatively large 
errors (20% for the 1-1 Pade). This indicates that the extrapolation procedure is probably 
far from optimal, and could be improved by using information on the asymptotic behavior 
of the large-?7 expansion for example. 

For finite A^, different Pade approximants give quite different answers, however it turns 
out that on the average the results are not worst than in the large N case. Results for 
the extrapolation of nip/Ajjg to m" = are presented in table || for different values of 
A^ and different Pade types, where we show the optimal values and the exact values from 
reference WM. In addition we illustrate in figure 5 a comparison of the a dependence for 



the different Pade types, for the lowest value of A^ = 2, the behavior for other values being 
very similar. 

We note the great stability as a function of a when a is varied in a wide range, in 
contrast with the dispersion of the curves in figure 0. We see that in some cases the 
optimal values are quite remarkably close to the exact ones, even for very low values of 



^^Note that the re-integration oi mp / Ajjg is unique, as the integration constant is fixed by the asymp- 
totic m" — > 00 behavior. 

^^The blank entries correspond to values of a where the coefficient Ci is positive. Ci is a complicated 
function of a, and there may be more than one range of values of a where it is negative, which happens 
for the first time for the 3-3 Pade. 
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1 0044 
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1 0080 


1 00Q3 

X . WWC/(J 




1 '^1 1 5 


1 0463 




n 7 


1 4S7Q 


1 0623 


Q526 


n 8 


1 51 (15 


1 0668 


Q5Q4 


n Q 


1 58Qfi 


1 06Q6 


Q602 

W. tJWWZj 


1 n 


X . U 1 uu 


1 070Q5 


Q582 

W. uOO^ 


1 1 

X . X 


1 771 1 

X.I 1 X X 


1 070Q6 


Q544 


X . ^ 


1 87fi8 

X . O 1 uo 


1 06Q8 

X . UUC/O 


Q4Q1 

W. tJHzC/ X 




2 (1025 


1 0676 


0.9422 


1.5 


2.2636 


1.0602 


0.9231 


opt. 


1.203 


1.071 


0.9602 








1.002 


error 


+20% 


+7% 


-4%, +0.2% 



Table 1: extrapolation of rrip/AjYs towards m" = using Pade approximants in the case 
= oo. 



A^. This is encouraging especially in view of the fact that the (only known) first two non- 
RG perturbative corrections were used. Actually however, the extremely good agreement 
for specific Pade types, for instance at A^ = 2 and A^ = 3, should be attributed to 
numerical occasion. What is certainly more significant is the average error over different 
values of A^, for a given Pade order and type. The overall best result for arbitrary A^ are 
obtained with the simple Pade (2,3) and the exponential Pade (2,2) which we note are 
both using expansion to the fifth order in l/rj in ( p.24| ). This is most probably due to 
the fact that, at the order where we are working, for a fixed A^ and a, expression (|3.14|) 
actually depends on five independent parameters: three perturbative terms, namely the 
"zero order" overall coefficient, Aii, and Ai2] P^us two independent combinations of RG 
coefficients, namely the coefficient of ln/2 in (|3.9|) and the power of /2 in ( p.l5|) . However 
it is not excluded that higher Pade orders would give even better results. Due to the rather 
complicated dependence of a given Pade upon the above parameters, we have not tried 
to systematically study this issue. It would certainly be of interest to further optimize 
the choice of the Pade approximants, but this would go beyond the scope of the present 
paper. At any rate, given the very different Pade approximant types confronted here, we 
can be confident that we have found a useful convergent variational scheme to compute 
non-perturbative quantities. 
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1 474QS iTl 4fi%'l 
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P2,3(m) 

opt. result (error) 


1.6206 (12.8%) 


1.3456 (9.2%) 


1.1917 (3.6%) 


1.1205 (1.1%) 


P2,3(«) 

opt. result (error) 


1.8749 (0.77%) 


1.4864 (0.3%) 


1.2654 (2.3%) 


1.1628 (2.6%) 



Table 2: Extrapolation of nip /Ajfg towards m" = using different Pade approximants 
for finite A^. The last two lines provide in addition a comparison of different RS for a 
same Pade approximant order: accordingly the results given in the last line were obtained 
using the alternative Ansatz form ( [4.14| ) derived in section 4. 



A particular feature of renormalizable theories might however spoil this optimism: it 
is the dependence of perturbative computations on the renormalization scheme (RS). As 
is well known, apart from the arbitrariness of the renormalization scale, parametrized 
with a, already at the second RG order there is more arbitrariness in the renormalization 
prescriptions. Typically, only the first two perturbative coefficients of the beta function 
and the first coefficient of the anomalous mass dimension are scheme independent, so 
that higher coefficients can be set to arbitrary values by perturbative redefinitions of the 
coupling constant and finite renormalizations. In our case, this translates into the fact that 
in principle the higher coefficients Aii for i larger than one could be set to arbitrary values 
with an ad hoc renormalization scheme. This problem is well known in perturbative QCD. 
The only assurance is that any reasonable renormalization scheme should give similar 
results at a given order, and that the accuracy of the results should improve as the order 
increases. What reasonable means exactly is unknown, but dimensional regularization 
with minimal subtraction is generally believed to be such a reasonable scheme. To test 
convergence and accuracy, we do not have higher order calculations at our disposal, but 
we can use the other method: do the same computation at the same order with a different 
renormalization scheme (RS) which can be considered a priori as reasonable, and compare 
the results. This is what we now do. 
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4 Alternative Scheme Ansatz 



The renormalization group invariant expression (|3.7|) has been reconciled in section 3 
with a variational mass expansion, using renormalization group invariance properties of 
appropriate combinations of bare parameters, mo, go like in equation (p.6|), and the resum- 
mation properties of the x series, resulting in the integral result ( p.8|) . However, once this 
construction is understood, and a finite, renormalized mass gap (|3.14|) emerges, one can 



alternatively try to construct an Ansatz starting directly from renormalized quantities]^. 
This will turn out to provide a more transparent intepretation of the different quantities 
involved in ( ^.8D -( P?T^ ), which also turns out to be more appropriate for a generalization 



to other theories, like QCD typically. 

Let us thus start from the renormalization group invariance properties of the renormal- 
ized mass and coupling constant, and examine what kind of non-perturbative information 
it may contain. Integrating the renormalization group equations for the running mass 
m{fi) in terms of the running coupling constant, g{fi), to two-loop renormalization group 
order exactlyQ, with the "fixed point" boundary condition 

one obtains after some algebra the expression 



RG 



1 + t9'm)f-' 



1 + 



Tl I 70 
■2bi ~^ 2bQ 



m{j2) f 

where / = g^{fi)/ g^{Mp'^) satisfies the recursive relation 



(4.1) 



f = l + 2bog\fi)ln^ 



1 + iro9'm) 



f 



(4.2) 



and the renormalization group coefficients 6j and 7j are given in Appendix B. Mp^ des- 
ignates the part of the mass which is entirely determined from renormalization group 
properties, i.e. it does not include the purely perturbative non-logarithmic finite parts. 
The latter are consistently included as follows: 



RG 



[I + mil — -r h m22- 



/ 



P 



(4.3) 



with mil and m22 simply given by the non-logarithmic parts of the perturbative expression 
for the renormalized mass, which is evidently obtained from equation ( p.l| ) by substituting 

^^This procedure will also be treated in more details in the QCD context [ p^ . 

^^i.e. keeping the exact dependence on g as given by the renormalization group but of course the 
renormalization group functions themselves being truncated at a given order. 
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niQ = Zmfn, Qq = (/i) ^ZgQ^. Using Z^, Zg at second order, one finds (in the MS scheme) 
Mr = m(^)(l-/(/.)^^^ln!^ (4.4) 



(iV-l/2)(iV-3/4) ^^,m(/i) 



^ (AT -l/2)(iV- 1/4)) ^^ m(/i) 
^ ^(0.737775 - — 



7r2 ^ 96 

from which we immediately obtain 



mil = , = ^ ,^^^ 0.737775 - ^) . (4.5) 

vr"' 96 



One can check that expansion ([4.3|) then reproduces the pertubative resuh (|4.4| ), and 
in fact generates correctly, by construction, the leading and next-to-leading logarithms, 
(yf^"ln"[m(/i)//i] and (7^" In""^ [m(/i)//i] respectively, to all orders. 

Next, we introduce the variational mass expansion in analogy with the procedure 
derived in section 2, by simply assuming the result to be given from performing the 
substitution 

m(/i) m{fi) u (4.6) 

everywhere in expressions ( [4.1| ), (^I2|), ( |4.3| ), and integrating the resulting expressions 
along an appropriate contour in the complex u plane, with the weight duu~^ e", similarly 
to e.g expression ( p.8|) . This leads to the final formula. 



2z7r Jc fnC + fl]B\ f, 
where we defined = f /{2bog'^{fi)) = [2bog'^{M^^) which satisfies 

/^(m) = \n[m"'u] -A\nf^-iB-C) \n[C + f^] . (4.8) 

In eqs. (|3) and (|]8D, A, B, C are defined as 

A = 2L = I1^ (49) 
26i 2{N-1)' ^ ' 



B 
C 



7^ _ 7^ ^ AT - 1/2 
26o 26i 2(A^ - 1) 
hi 1 
26g 2(A^ - 1) ■ 



and 



= U 10) 

A^2-^[26o^;^(/.)]^+^[l + f^^?2(^)]-i.- ^ • ) 
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In equation ( |4.7|) , we also included the RS dependence through the parameter a, in a way 
similar to (|3.14| ). Its explicit dependence, dictated by the renormahzation group, reads: 



M\{a) 



(N - 1/2) 



(N-r 

(N - 1/2) 



In a , 



(4.11) 



(iV- 



(0.737775 



-l/2)(iV 



1/4) 



(N-l) 



96' + 
In a . 



(A^- 1/2) (A^- 3/4) 
(AT - 1)2 



In a 



It is instructive to examine more closely the differences between the two forms ( p.l4|) and 
( [4.7| ): at first sight they look very different, as not only the functional form of Mp, the 
power coefficients, and the perturbative coefficients A4i, A4[ are different, but the form 
of the recursive functions, /2 in equation (|3.15|) and in equation (^ 



are also different. 



Actually, we shall see next that this is in fact nothing but a RS difference, although in 
a not too conventional form. Since, at the second RG order, the RG properties allows such 
RS changes, this proves a posteriori that the construction when substitution ( [4. 61) is made 
in the relevant renormalized expressions, followed by the appropriate contour integration, 
gives an equally acceptable Ansatz. In particular it corresponds, up to RS change, to an 
expression resumming infinite orders of the x series expansion, as is originally aimed. 

To make the precise connection between the two different RS choices, note first that 
the definition of the variational parameter m'" in ( 4.10| ) and m" in equation ( |3.14 ) indeed 
differs (except for iV — >^ oo): 



m 



exp 



-1 



2(iV- 1) 



m 



(4.12) 



The latter relation, by also comparing Q the form of ( [4.8| ) with /2 in ( p.l5| ) allows to 
express perturbatively the relation between the two schemes, for and /2 sufficiently 
large: 

1 (3iV + 1/2) 
J2 ^ 



1 



1)^ /I ■ " VI 



(4.13) 



2(iV-l)/2 A{N 

One can easily check, by expanding in equation (|4.7| ) in powers of I//2 using ( [4.13|) , that 
one recovers, up to higher order 0(1//|) terms, the previous form ( |3.14| ) Consequently 
the two schemes are perturbatively equivalent, as they should. For f2 —>■ (equivalently 
— i> 0), i.e. m" —>■ (equivalently m"' —>■ 0) the two different forms have however no a 
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From ( [4.12D one may think at first that this RS change is simply given as a specific choice of a. This 



is not so, since the different form of in (4.8) with respect to (3.15) is also essential, especially for 
m'" 0, and cannot be obtained from the alternative form /2 by simply changing a. 

But with 62 = 72 = 0, because in the alternative Ansatz these third order RG coefficients do not 
appear, while in the previous Ansatz those are needed at intermediate stage for consistency, see Appendix 
B for more details. 
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priori reason to give the same results, due to a different behavior of /2 and close to the 
origin. 



To proceed, we perform the same type of Pade approximants, as described in section 
3, with the new form of the mass gap ( [4.7|) . Actually formula ( [4 .71) as it stands is not 
directly convenient, due to the presence of an extra branch cut starting at /2 = —bi/ (26q) 
i.e. on the positive real axis {bi < 0), therefore preventing a continuation down to m'" —>■ 
like in the calculations of section 3. One may possibly infer that 6i < is an artifact 
of the low orders of perturbation theory {i.e. the Gross-Neveu model is asymptotically 



free to all orders in the large- limit [TE\). But in any event the analyticity structure is 
not uniquely determined, since it clearly depends on the scheme and precise form of the 
defining relation for /2, as illustrated here from the two different expressions (|3.15|) and 
([4 .81). It is in fact possible to recover a structure of singularities similar to the previous 
case (while still using the different information provided from the alternative scheme) by 
simply expanding (|4.71 ) for large /a (namely in the perturbative range above the extra cut 
at /2 = 1/{2{N — 1)), but keeping the same definition of the variational parameter m'" in 
equation ( |4.10| )["]. This replaces expansion ( |4.7| ) by 

mp ^ 1 m'" a r due"" r M"i{a) M"2(a) } , , ^ 

where /2 verifies the very same relation as in equation (|3.15|) (but now with m" m'"), 
and the cut again starts at /2 = towards the real negative axis. The perturbative 
coefficients are of course different: 



-M"i(l) 



A^-1/2 
4(iV- 1)2 



(IQN^ - 14A^ - 3) 
^"2(1) = A^-(l)-(iV-l/2) ^ 64(iV-l)^ • ^^-^^^ 

One can now directly apply the Pade techniques as described in section 3, to the form 
( [4.141) . The most interesting results are illustrated with the Pade (2,3) case in the last 
line of table 2, and are substantially better than the corresponding order results in the 
previous scheme. In particular the error is much more stable when varying N. Optimal 
values in this alternative scheme for Pade approximants of lower orders are only slightly 
different from those of the original scheme as given in table 2. For larger A^, the two 
schemes give more and more similar results, as expected since for N 00 they are 
strictly equivalent: m'" m", and all perturbative coefficients A4i{l) go to 0. These 
different facts provide a good check of the stability of our results with respect to changes 
in the renormalization scheme. Indeed it strongly indicates that it should also be possible 

^''Of course the larger N is, the closer one can safely approach the origin for (equivalently m'"), as 
expected. 
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to optimize with respect to the RS scheme (in addition to the optimization with respect 
to a). The variation of the scheme through the parameter a (whose dependence does 
differ in the alternative Ansatz at the next-to-leading order, (lna)//|), provides an extra 
illustration of the consistency of our results. Besides, these two different renormalization 
schemes are important to clarify the variational expansion procedure when applied either 
to bare or renormalized parameters. To illustrate a more general RS dependence one could 
further exploit the RS dependence of the renormalization group coefficients themselves, 
e.g. starting with 71. Although we have not tried explicitely this in the present case, we 
see no compelling reasons to study a more general case within the framework of the GN 
model, where we know the exact results anyway. The preceding facts should already be 
considered a strong indication that the variational mass method does lead to non-trivial 
results, at least in appropriate renormalization schemes. 

The vacuum energy density can be treated along the same lines. However, in that 
case, the change of sign of the energy when m" approaches complicates matters in a 
way which we do not report here, as it is a separate issue. 



5 Conclusion 

In this paper, we have shown how to introduce a variational procedure in an asymp- 
totically free field theory, in a way which is compatible with the renormalization group: 
we have found how to rescale the variational parameter in order to take account of all 
infinities (at least in the framework of dimensional regularization) so that the physics of 
the variational approximation is smooth when the cut-off is removed (in our case, when 
the dimension is varied around its critical value). We have further shown that the nu- 
merical results obtained with an extra very simple extrapolation procedure are in good 
to remarkable agreement with exact values in the case of the GN model. This establishes 
the theoretical possibility and potential usefulness of perturbative calculations for com- 
puting non-perturbative quantities. This framework is in principle directly applicable to 
the QCD case, apart form possible complications due to the a priori different analyticity 
structure of m of the relevant QCD expressions. In a subsequent paper [Q, using a bare 



chiral symmetry breaking fermion mass as variational parameter, we shall obtain values 
for some parameters of spontaneous chiral symmetry breaking, mq (the part of the con- 
stituent quark mass due to chiral symmetry breaking), the pion decay constant and 
the condensate < ipip >. 
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A A contour integral resumming the x dependence 

We shall explain here in some details how to resum the x series generated from the 
perturbative expansion of ( |2.1| ). Consider the one-loop RG invariant (bare) expression for 
the mass mi;', as given in eq. (|2.3|): 



rriF = z (A.l) 

where the substitution 

mo mo(l - a;); gl glx, (A.2) 

provides a new quantity m^(x). To pick up the order term in mp{x) = I]^o S-^"^ 

(having in mind that we are actually interested in the limit x ^ 1), a. convenient trick is 
by contour integration: 

mf[x 1) E «fc = 7^ / dx{- + ■■■ + mF{x). (A.3) 



Now performing the sum in ( |A.3|) exhibits a (1 — x) ^ factor which cancels the (1 



X] 



from (|A.2| ). This results in the expression 

mf{x 1) ^ — (f dxx-'-'+^^ rno[fo{x)r^ , (A.4) 
2z7r J 

where the contour is counterclockwise around the origin, and for convenience we defined 
the (recursive) function 

foix) = 1 - Anbo X glT m^(l - xY (/o)"^^ , (A.5) 
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In (A.4) there appeared in fact a factor of 1 — x (9+^), from which only the last term contributes to 



the integral due to the analyticity of fo{x) defined in (A.5) 
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directly dictated from eq. ( |A.1| ). /o(x) has (evidently) a power series expansion in x, but 
admits also an expansion in (1 — x), as noted by inverting its defining relation (|A.5|) . This 
implies, in particular, that x = 1 is an (isolated) pole of mi;'. 

Provided that no extra singularities lie in the way, one may distort the integration contour 
in (|A.4|) to go around the cut lying along the real positive axis and starting at x = 1. 
Actually one can go a step further and reach the q oo limit: after distorsion of the 
contour only the vicinity of x = 1 survives for g ^ oo, that one can analyse by changing 
variable to 



V 



l-x=- (A.6) 

q 

and rescaling mo by introducing mo = mgg (keeping mg fixed as q goes to infinity). One 
finds in place of ( |A.4|) 

.4"to-oo)^^/*e"-^ (A,7) 



V 



fo{v) 



where now /o(^^) = 1 — 47r6o 5'0^r(mof )''(/o) "^'o . The crucial point in ( |A.7|) is, that once 



performing renormalization via mo = mZm, Qq = fi ^ZgQ^, nip is finite to all orders: 

1 / m 

mp = f dve — (A. 8) 

2z7r J j2b^ 



where the renormalized function / = Zgfo = 1 + 26ofi'^ ln[(mf //i) / 



-(70/260) 



We have thus recovered finite quantities with a non-trivial x expansion. In the latter 
derivation we only included the one-loop RG dependence, which is the exact result in 
the large limit only. For arbitrary one should include in the derivation the non- 
logarithmic perturbative terms, present e.g at the two loop order in (|4.4| ). This can be 
done without affecting the contour integration properties, except that the resulting ex- 
pression of rriF has a more complicated structure around f ^ 0, but can be systematically 
expanded around the origin in the way discussed in section 2. Generalization of the pre- 
vious construction to the next RG order is straightforward, since the recursive function / 
has a very similar form as in ( |A.5| ) above, where only the power coefficients are changed. 



B Bare RG Ansatz 

We give here some useful expressions needed for the construction of the generalized Ansatz 
using third order RG functions, as introduced in section 2.3, and leading to our main 
results (p. 81), ( |3.9| ) and ( |3.14| ) for the mass gap. Derivation of the vacuum energy ( p.l7| ) is 



very similar, apart from the slight complication due to the subtracted terms as discussed 
at the end of section 2. 
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We start from expression ( p.7|) , where the constants A, B and D are first to be deter- 
mined by requiring equation (|3.71 ) to be finite as e — >■ 0, when expressed in terms of m 
and with the help of RG mass and couphng counterterms at a given order (and in a 
specific renormahzation scheme). The coefficient X, p, . . . are then obtained by matching 
the perturbative expansion of (|3.7| ) to the perturbative (bare) expression ioi mp, eq. ( pTII) . 

Defining the renormahzation group coefficients governing the couphng constant and 
mass evolution, respectively, as 

P{g) ^p^='-g-bog'-hg'-hg'~... (B.l) 



and 



^-^^ = lo9' + li9' + l2g' + ... (B.2) 
m dfi 



the coefficients bi and 7^, known in the MS scheme up to three loop order in the 
Gross-Neveu model, are expressed as 

'''''^ N-l/2 AT -1/2 (iV-l/2)(4iV-3) 

From these one obtains explicit expressions of the counterterms, Zg and Z^, as expansions 
both in g^ and e to the required order. We need a counterterm for in ( |3.7| ) as well, 
that we define according to jFg = ZpJ^, where Zp can be determined consistently together 
with A, B and D for given Z^ and Zg. In terms of Zg and Z^, the finiteness of eq. (|3.7|) 
requirement explicitely reads 

B lnz + lnZm= finite; (B.5) 
(1 + eA) Inz + e In Z^ + In Zg = finite; 
D = l + eA- 



where we exhibited the singular part of Zp as 



In 2; = In 



9' {N - l)g' 



1 + ^- 



[B.6) 
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Prom inspection of the formal expansion of ( PTT] ) it turns out that actually A is needed 
to 0{e) and B to 0(£^), to fix unambiguously the finite coefficients A and p in (|3.3| ) by 
matching with (|3.1| ). Accordingly from ( [B.5|) one needs to know Zm, Zg to order 0(1/ e^): 
i.e the third order RG coefficients 62? 72 in (|B.3| ), ( |B.4| ) should be included for consistency. 



The required expansion in powers of 1/e of InZg, In Zm and \nz is completely deter- 
mined by the expansion in e of the perturbative zero, gp, of the beta function, which is 
easily derived as 
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In this way we obtain after straightforward algebra, 

N N -1/2 , , 

A^-1/2 , 3Ar-l/2 _ N{N -1/2) 3 



and0 



^ = lV^ + 4(iV^^+ 8(iV-l)3 ' +^^^ 



477 - 1 ^ ' 



+ 



iV- 1/2 
B7r(iV- 1) 



N 



+ 3(7i5-ln47r) 



(A^ - 1/2)(36A^2 _ g2]v + 17) 



2 



647r2(iV - 1) 
' ' -0.725551... + — +0 e 



27r2 ' 48' 

which directly leads to expressions (|3.9| ) for e — 0. Let us finally remark that ( |B.5| ) leaves 
over non-zero finite perturbative terms, which are however completely determined order 
by order. This is inherent to our scheme here where the counterterm for T could not be 
put into a form where only the singularities appear (like is the case, by definition, for 
InZjn and InZg in the MS scheme). Explicitely after using ( p.8| ) we obtain 

(A^ - 1) 47r (A^ - 1) 167r2 ^ ^ ' ^ ' 



which obviously gives perturbative corrections when re-expanding expression (|3.8| ), and 
which are indeed necessary to recover consistently the perturbative expression for the 
renormalized mass ([4.4|) 0. 



^^Note that the order e term in A is required to determine p. 

^°The third order renormalization group dependence, via &2 and 72 in (B.8), ( |B.g| ) and ( B.10| ) actuaUy 
cancels in the (two-loop) perturbative renormalized mass expression (4.4). 
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Figure Captions. 



Figure 1. Mass gap graphs at order 2. 
Figure 2. Vacuum energy graphs at order 2. 

Figure 3. Mass gap ratio m.F/A.jjg versus variational parameter m", for values of the 
scaling parameter a ranging from a — 1.1 (lowest curve) to a = 1.2 (highest curve). The 
horizontal dashed curve represents the exact value. We have taken N — 3 ior the number 
of flavors. 

Figure 4. Vacuum energy Eq/Kjj^^ versus m", for values of a ranging from a = 0.9 (lowest 
curve) to a = 1.1 (highest curve). The horizontal dashed curve gives the exact value. The 
number of flavors is = 3 . 

Figure 5. Comparison of different Fade approximant types and orders, versus exact 
mp/^-Ms^ for N — 2. The curves showing extrema at 2.27, 1.62, 1.76, 1.93 and 1.87, 
have been obtained from the pi^i{u), P2,3{u), P2,2{u) x e^'^^^^\ Pi,2(it) and p^^^iu) Fade 
approximants, respectively. 
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Figure 1: Mass gap graphs at order 2. 
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Figure 5: Comparison of different Fade approximant types and orders, versus exact 
mp/AMS, for = 2. 
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